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Abstract. The upper Asymptotic Giant Branch (AGB) is populated with oxygen rich and carbon rich Long 
Period Variables (LPVs). These stars are essential contributors to the near-IR light of intermediate age stellar 
populations. Individual observed spectra of LPVs are so diverse that they cannot be used directly in the synthesis 
of galaxy spectra. In this paper, the library of individual spectra of Langon & Wood (2000) is used to construct 
averages that can be incorporated conveniently in population synthesis work. The connection between such spectra 
and stellar evolution tracks is discussed. In order to select a sorting criterion and to define averaging bins for the 
LPV spectra, correlations between their spectrophotometric properties are reexamined. While optical properties 
and broad baseline colours such as (I-K) are well correlated, a large dispersion is observed when these indices are 
plotted against near-IR ones. This is partly due to the intrinsic width of the upper AGB, which is illustrated by 
locating each of the multiple observations of individual LPVs on the HR diagram. It is argued that broad baseline 
colour-temperatures are the most sensible sorting criteria. The properties of the resulting sequence of average 
spectra indeed vary regularly. We further address: (i) the bolometric corrections and temperature scales needed to 
associate a spectrum with a given point on a theoretical stellar evolution track (or isochrone) , (ii) the simplifying 
assumptions that will be implicitely made when using the average spectra, (iii) potential biases in the sample of 
Langon & Wood and their effects, (iv) the small contribution of LPVs to the interstellar hydrogen emission lines 
in galaxies. It is emphasized that an a posteriori calibration of the effective temperature scale remains necessary, 
until consistent models for the evolution, the pulsation and the spectral appearance of LPVs become available. 
We suggest a recipe for the use of the average spectra at various metallicities. 

Key words. Stars: AGB and post- AGB - Stars: late-type - Stars: variable: general - Hertzsprung-Russell diagram 
- Infrared: stars - Galaxies: stellar content 



1. Introduction 

The prediction of the near-infrared (near-IR) spectra of 
star clusters and galaxies is subordonate to the existence 
of stellar spectral libraries with a complete coverage of the 
evolved stages of stellar evolution. Stellar libraries are also 
necessary in the construction of near-IR colour-magnitude 
or two-colour diagrams, which may include narrow-band 
filter indices or absorption line measurements. A short- 
coming of all libraries previously us ed in population syn- 
thesis work (e.g. Kleinmann & Hall |l986 Terndrup et al. 



1991] , Langon & Rocca-Volmerange |199J, Pickles |l998| ) is 



the poor coverage of the upper Asymptotic Giant Branch 
(AGB). The upper AGB hosts oxygen-rich and carbon- 
rich Long Period Variables (LPVs)[], whose specific spec- 



Send offprint requests to: A. Langon 

1 In agreement with common usage, LPVs include both 
semi-regular and Mira-type variables; Miras are LPVs with 
large (opti cal) amplitud es, i.e . SV>2. 5 magnitu des (K holopov 
et al. [1984 L1 °y d Evans Ejjgl Hughes & Wood |l99C|). 



tral features carry the potential of revealing the pres- 
ence of intermediate age populations (Langon et al. 1999). 
Indeed, upper AGB stars contribute of the order of 50 % of 
the K band light over a range of stellar population ages be- 
tween ~ 10 s and ~ 10 9 years (Persson et al. 1983, Rcnzini 



& Buzzoni 1986, Frogel et al. 199C| , Ferraro et al 



1995 



Bressan et al. 1998, Girardi & Bertelli 1998, Langon 



Maraston |1998[ Mouhcine & La ngon |2002j ) 



1998 



Recently, Langon & Wood (|2000j hereafter LW2000) 
published a library of stellar spectra of luminous cool stars 
that includes a large sample of instantaneous observations 
of LPVs. The data span the wavelength range where most 
of the light of cool stars is emitted (0.5 to 2.5 jim) with a 
spectral resolving power of ~ 1100 longward of 1 /^m, and, 
with about 100 spectra, provide the most complete data 
set available to date. 

The range of spectrophotometric properties observed 
in LPVs is large (LW2000, Alvarez et al. |2000[ ). Even the 
classification of the empirical spectra has many caveats. 
Sorting algorithms based on optical or on near-IR criteria, 
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on broad band colours or on spectroscopic properties, on 
mean or on instantaneous data, often give different results. 
The most conspicuous features, such as the H2O vapour 
absorption bands around 1.4 /zm and f.9 fim (in oxygen 
rich LPVs), depend on details of the atmospheric structure 
that change from one pulsation cycle to the next. Such 
changes will be smoothed out in the integrated light of 
(large) stellar populations. 

The purpose of the present paper is to provide a library 
of mean stellar spectra of luminous cool variable stars, 
suitable for direct use in combination with a population 
synthesis model. The library complements more readily 
available data for static red giants or supergiants. Both 
oxygen rich and carbon rich AGB stars are considered. 
We attempt to thoroughly discuss the impact of techni- 
cal/physical choices that have to be made at various steps 
in the preparation and the use of the averaged spectra. 

The construction of a library of mean spectra requires 
the preliminary construction of suitable averaging bins. 
In principle, one would like to compute energy weighted 
mean spectra for individual stars, based on spectra taken 
at various phases over several pulsation cycles and on the 
light curve. However, this requires an enormous amount 
of spectroscopic and photometric data. Even the LW2000 
data are not sufficient. This paper suggests that colour 
temperature is the best single parameter to order LPV 
spectra. The sample of useful input spectra is briefly de- 
scribed in Sect.|J where the main selection criteria (and 
thus potential selection biases) are also recalled. Sequences 
of average spectra are presented in Sect.||. The statistical 
properties of the data that justify our selected sorting cri- 
teria are also presented in that section. 

The discussion of the uncertainties in population syn- 
thesis predictions inherent to the chosen averaging proce- 
dure must be based on the understanding of two processes 
that tend to disperse the stars of even a coeval population 
over quite a broad area of the Hertzsprung-Russell (HR) 
diagram: the TP-AGB thermal pulse cycles and the LPV 
pulsation cycles. The theoretical stellar evolution tracks 
used in population synthesis calculations provide the time 
evolution of the position of a star in the HR diagram. The 
thermal pulses are now commonly accounted for in the 
tracks, but the effects of pulsation, on timescales of hun- 
dreds or thousands of days, arc not included. One may 
think of these tracks as providing the evolution of the 
static parent stars of LPVs. Unfortunately, building non- 
linear pulsation models and theoretical spectra for LPVs 
is extremely complex, and current models cannot yet pro- 
vide the fundamental parameters of the parent star on the 
basis of one or several instantaneous empirical spectra or 
colours. 

How, then, should one connect the averaged spectra 
of the new library with locations along the evolutionary 
tracks? Practical aspects of this question are addressed in 
Sect.|[ The subsequent discussion addresses fundamen- 
tal difficulties. In particular, it emphasizes the need for a 
posteriori calibration of the relation between colour tem- 



peratures and effective temperatures (i.e. the tempera- 
ture scale of the spectra) . Potential selection biases in the 
LW2000 data, and their effect, are also discussed in Sect. [| 
Recipes for extensions of the population synthesis calcu- 
lations to other metallicities than quasi-solar are provided 
in Sect.|^. A concluding summary is given in Sect.^. 

The averaged spectra are available in digital form 
through CDS0. 



2. The sample 

The library of LW2000 contains more than 100 merged 
optical/near-IR spectra. Of these, about one quarter be- 
long to giants of the first giant branch, red supergiants, 
Galactic Bulge LPVs or luminous red stars of the Large 
and Small Magellanic Clouds. The remaining spectra are 
those of pulsating giants in the field of the Milky Way. 
These are the data we focus on. 

The sample of LW2000 was selected to provide a rela- 
tively uniform coverage of the period-amplitude plane of 
LPVs (their Fig. 1). The O-rich LPVs with V band am- 
plitudes (<5V) smaller than 0.7 magnitudes in the sample 
show no obvious difference with static giants. We exclude 
these objects from the subsample used here. The resulting 
set contains 63 O-rich spectra of LPVs with SV > 0.8 mag. 

The sample of C-rich stars contains only 6 objects, 
of which one is an S/C star (3 of 21 spectra) and one 
is a far-IR source (R Lep = IRAS 04573-1452; 3 spectra). 
Although their amplitudes are larger than 1.5 mag in V, 
the spectra show relatively little variation with phase. 
The largest amplitude amounts 5V~6mag for the far-IR 
source. The next largest amplitude is <5V~3.7. The spec- 
trophotometric properties of the C stars in the sample are 
much less dispersed than those of the O-rich LPVs. 



3. Average spectra for LPVs 

In principle, the ideal spectral library for population syn- 
thesis purposes should consist of a large collection of light- 
curve weighted spectra of individual stars observed pho- 
tometrically and spectroscopically over several pulsation 
cycles. There are two main reasons why this is not prac- 
tical. First, the amount of data required is enormous. In 
the sample of LW2000, only a few light curves have been 
sampled at more than 4 points. The actual light curves 
and phases are generally not accurately known. No more 
complete empirical spectral library for LPVs is available 
in the literature. Second, uncertainties in current models 
make it impossible to determine the physical nature of one 
particular star (initial mass, current evolutionary status), 
even when many spectra are available for it. Attempts 
to follow the "ideal" route fail at a very basic level: the 
mean spectra of individual stars cannot be sorted into a 
sequence that behaves regularly enough for practical use. 



2 Centre de Donnees Astrophysiqu es de Strasbourg, 
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Along the TP- AGB, the effective temperature globally 
decreases while the luminosity increases. In this section, 
we show that the use of effective temperature indicators as 
basic sorting criteria produces a sequence of O-rich spec- 
tra with regularly evolving properties. The tightness of 
this empirical sequence tells us that surface temperature 
is indeed a good parameter for the first order classification 
of the spectra. Other sorting parameters are discussed but 
are found to be less appropriate. Using the suggested tem- 
perature sequence is, in our opinion, the most sensible way 
to include O-rich LPVs in evolutionary synthesis models 
for galaxies. 

In the case of C-rich stars, both a temperature se- 
quence and a sequence of C/O abundance ratio are con- 
sidered. 

3.1. Oxygen-rich LPVs 
3.1.1. Temperature sequence 

The spectra of oxygen-rich LPVs selected in Sect.|| cover 
a range of 8 magnitudes in (V-K) and 4 magnitudes in (V- 
I) or (I-K)[J These long baseline colours are expe cted to 



corre late with effective temperature (Bessell et al. 1989 b 
1998| ), although the absolute calibration of the relation 



is poorly known for variable stars, and is metallicity- 
dependent (cf. Sect.||). 

Figure [j] confirms that tight correlations exist between 
various commonly used photometric temperature indica- 
tors. The small dispersion in panels (b) and (c), which 
relate (V-K) and (R-K) to (I-K), shows that the three 
long baseline broad band colours contain equivalent infor- 
mation. In the following, (I-K) is adopted as the primary 
temperature index. The choice between (I-K) and (R-K) is 
arbitrary, while (V-K) was rejected because the passbands 
of standard V filters extend significantly over the edge of 
the LW2000 data, requiring extrapolation of the data or 
(as was done here) a truncation of the filter passband. 

The VO index of panel (f) measures the depth of the 
VO A-X (Av — 0) band at 1.05 ^m, as done by Bessell 
et al. ( [1989 b| ) and Alvarez et al. ( |2000D . The absence of 
the molecular band in the warmer LPV spectra restricts 
its use as a temperature indicator. S1/3 and S2/3 in pan- 
els (d) and (e) are the optical narrow band filter indices 
defined by Fluks et al. ( |1994[ f|. For static giants, they are 
good indicators of the spectral type, although the molec- 
ular bands that determine S1/3 saturate for the coolest 
stars (spectral types later than M6). The relative disper- 
sion in the figures involving narrow bands is due to the 



3 Unless otherwise stated, broad band colours are measured 
on the reduced spectra using the transmission curves of Bessell 



(L99C) for V, R and I (Cousins system), and of Bessell & Brett 
(L988) for J, H and K (Johnson-Glass system). Zero colours 
are adopted for Vega, based on the Vega model spectrum also 
used by Bessell & Brett (we thank M. Bessell for providing the 
latter) . 

4 Note that integrated flux ratios are used in that definition, 
instead of the usual flux density ratios. 



sensitivity of individual molecular bands to the structure 
of the outer atmosphere, which is affected by details of 
pulsat ion. Th is sensitivity has been discussed by Bessell 
et al. 1989 b on the basis of Mira models and is found 
empirically in measured variations of the bands relative 
to each other or to the colo urs, fr om cycle to cycl e and 
with phase (Spinrad & Wing 196E, Alvarez & Plez, 1998 , 
LW2000). It strengthens the preference given to (I-K) as 
a temperature indicator. 
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Fig. 1. Temperature sensitive indices as a function of the 
bin number based on (I-K). The indices of the O-rich LPVs 
of LW2000 are shown as crosses; the solid lines connect 
those of the averaged spectra. All colours are in magni- 
tudes, with a value of for Vega. In (b), V is measured 
through a standard filter truncated at 5100 A(cf. Alvarez 
et al. |200C|). 



(I-K) was used to sort the 63 individual spectra into 9 
bins of 7 spectra each, as listed in Table]!] and illustrated 
in Fig. 0(a). As a result of this choice, bin sizes in terms of 
(I-K) are not constant. No binning choice guarantees con- 
stant bin sizes in all colours, and choosing constant bin 
sizes in (I-K) rather than a constant number of spectra 
per bin would have given (I-K) more importance than it 
deserves compared to the other colours of Fig.[l]. Within 
each bin, the spectra were normalised to a common mean 
flux over the whole wavelength range of the data, and were 
then averaged without further weighting. The resulting se- 
quence of average spectra is displayed in Fig.^. The maxi- 
mum of the energy distribution is seen to move from about 
7500 A for the hottest spectrum to about 1.3 /zm for the 
coolest. Molecular bands grow more numerous and deeper 
with decreasing colour temperature. The corresponding 
photometric index measurements are located on the solid 
lines in Fig.]!]. 



4 



A. Langon & M. Mouhcine: Average Spectra for Upper AGB Stars 



Table 1. Contents of the 0-rich temperature bins 



Bin 


Name 


Date 


Type 


Period 


5V 


(I-K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


TCen 


Jul. 96 


Sa 


90.4 


3.5 


1.86 


1 


TCen 


Jan. 96 


Sa 


90.4 


3.5 


1.90 


1 


BD Hya 


Jul. 96 


Sa 


in A 


1.8 


1.90 


1 


SCar 


Jun. 95 


M 


149.5 


5.4 


1.90 


1 


SCar 


Dec. 95 


M 


149.5 


5.4 


1.98 


1 


TCen 


Mar. 96 


Sa 


90.4 


3.5 


2.01 


1 


RS Hor 


Mar. 96 


M 


202.9 


5.2 


2.19 


2 


SLib 


May 96 


M 


192.9 


5.5 


2.26 


2 


SLib 


Jun. 95 


M 


192.9 


5.5 


2.29 


2 


SCar 


May 96 


M 


149.5 


5.4 


2.32 


2 


TCen 


May. 96 


Sa 


90.4 


3.5 


2.36 


2 


UCrt 


Jan. 96 


M 


169. 


4. 


2.45 


2 


BD Hya 


May 96 


Sa 


117.4 


1.8 


2.53 


2 


BD Hya 


Jan. 96 


Sa 


117.4 


1.8 


2.61 


3 


SCar 


Jan. 96 


M 


149.5 


5.4 


2.71 


3 


RPhe 


May 96 


M 


269.3 


6.9 


2.72 


3 


SCar 


Jul. 96 


M 


149.5 


5.4 


2.73 


3 


SPhe 


Dec. 95 


SI) 


141. 


2. 


2.75 


3 


UZHya 


Jan. 96 


M 


261. 


5.7 


2.78 


3 


UZHya 


Mar. 96 


M 


261. 


5.7 


2.82 


3 


RPhe 


Jul. 96 


M 


269.3 


6.9 


2.83 


4 


UCrt 


May 96 


M 


169. 


4. 


2.93 


4 


KV Car 


May 96 


Sb 


150. 


0.8 


2.94 


4 


RS Hor 


Jan. 96 


M 


202.9 


5.2 


2.98 


4 


RSHya 


Jun. 95 


M 


338.6 


5.2 


3.19 


4 


RSHya 


May 96 


M 


338.6 


5.2 


3.22 


4 


SPhe 


Jan. 96 


SI) 


141. 


2. 


3.24 


4 


SYVel 


Mar. 96 


Sb 


63. 


1.3 


3.27 


5 


KV Car 


Jan. 96 


Sb 


150. 


0.8 


3.31 


5 


SCar 


Mar. 96 


M 


149.5 


5.4 


3.38 


5 


SYVel 


Jan. 96 


Sb 


63. 


1.3 


3.45 



Bin 


Name 


Date 


Type 


Period 


5V 


(I-K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


5 


SPhe 


Jul. 96 


Sb 


141. 


2. 


3.49 


5 


SPhe 


May 96 


Sb 


141. 


2. 


3.51 


5 


RYCra 


Jun. 95 


M 


195. 


1.9 


3.51 


5 


RSHya 


Jul. 96 


M 


338.6 


5.2 


3.54 


6 


RS Lib 


May 96 


M 


217.6 


6. 


3.63 


6 


RSHya 


Apr. 96 


M 


338.6 


5.2 


3.70 


6 


RZ Car 


May. 96 


M 


272.8 


6.2 


3.73 


6 


SVTel 


Jun. 95 


M 


225.5 


3. 


3.76 


6 


RPhe 


Dec. 95 


M 


269.3 


6.9 


3.77 


6 


RS Hor 


Dec. 95 


M 


202.9 


5.2 


3.97 


6 


RCha 


Dec. 95 


M 


334.6 


6.7 


4.00 


7 


RSHya 


Mar. 96 


M 


338.6 


5.2 


4.07 


7 


XMen 


Jul. 96 


M 


380. 


3.6 


4.07 


7 


RPhe 


Jan. 96 


M 


269.3 


6.9 


4.23 


7 


RCha 


May 96 


M 


334.6 


6.7 


4.34 


7 


WWSco 


May 96 


M 


431. 


4.1 


4.36 


7 


RZ Car 


Jun. 95 


M 


272.8 


6.2 


4.49 


7 


RZCar 


Jul. 96 


M 


272.8 


6.2 


4.49 


8 


SVLib 


Jul. 96 


M 


402.7 


1.3 


4.56 


8 


SVLib 


Jun. 95 


M 


402.7 


1.3 


4.65 


8 


SVLib 


May 96 


M 


402.7 


1.3 


4.70 


8 


RCha 


Jan. 96 


M 


334.6 


6.7 


4.86 


8 


CM Car 


Jan. 96 


M 


335. 


2.5 


5.05 


8 


RSHya 


Jan. 96 


M 


338.6 


5.2 


5.24 


8 


XMen 


Dec. 95 


M 


380. 


3.6 


5.32 


9 


RZCar 


Mar. 96 


M 


272.8 


6.2 


5.38 


9 


XMen 


Mar. 96 


M 


380. 


3.6 


5.46 


9 


XMen 


Jan. 96 


M 


380. 


3.6 


5.55 


9 


RZCar 


Jan. 96 


M 


272.8 


6.2 


5.57 


9 


WWSco 


Jun. 95 


M 


431. 


4.1 


5.61 


9 


RSLib 


Jun. 95 


M 


217.6 


6. 


5.79 


9 


RCha 


Mar. 96 


M 


334.6 


6.7 


5.95 



3.1.2. Near-IR properties 

The spectral features of highest relevance to extragalactic 
studies are located longward of 1 /xm, i.e. at wavelengths 
where the relative contribution of the AGB stars to the in- 
tegrated light can be significant. The VO band of Fig. | (f) 
is one example. The behaviour of the near-IR broad band 
colours (J-K), (H-K) and (J-H) and of the CO and H 2 
absorption bands is shown in Fig.||. All these indices are 
based on the Hubble Space Telescope NICMOS filter pass- 
bandsf]. They are given in magnitudes, with the assump- 
tion that all colours are for Vega. The solid lines refer 
to the properti es of our prefered sequence of averages, de- 
scribed in Sect |3.1.l| . The dashed and dot-dashed lines are 
obtained respectively with 7 bins of 9 spectra each, and 
with 10 bins of 6 or 7 spectra each. 

Fig. U shows that the near-IR broad band colours ( J- 
K), (H-K) and (J-H) correlate with the (I-K) colour- 
temperature, but with a large dispersion and a relatively 

5 Many molecular index definitions are found in the liter- 
ature. We recommend that users measure the index of their 
choice on the actual spectra instead of attempting a conver- 
sion from the plots in this section. 



shallow slope. // (I-K) really is a good effective temper- 
ature estimator for O-rich LPVs, then it is clear that 
the near-IR colours are poor ones. Bins based on near- 
IR colours would mix spectra with very diverse general 
energy distributions, optical spectra and spectral types. 

The observed decrease in the near-IR colour indices be- 
tween (I-K)=3.3 and (I-K)=4.5 (bins 5 to 7) is the result 
of the competition between the influences of temperature 
and of molecular absorption. This decrease is statistically 
significant in our sample, in particular in (J-H): the size 
of the step is about 3 times the estimated r.m.s. error on 
the mean colour (r.m.s. dispersion of the data divided by 
v7j as we have 7 spectra per bin). However, it is unclear 
whether it would be seen in larger samples. In addition, 
the step is relatively small when compared to the whole 
range of the data, and it is comparable to typical observa- 
tional uncertainties in near-IR photometric observations. 
Our summary of these data is that the mean (J-K), (H- 
K) and (J-H) of O-rich LPVs are approximately constant 
between (I-K)~3 and (I-K)~5 (i.e. bins ~3 to ~8). 

Spectra with similar energy distributions display an 
impressive variety of near-IR molecular absorption spec- 
tra, as illustrated in Fig.0. Again, this shows how sen- 
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Fig. 2. The sequence of O-rich LPV spectra, obtained af- 
ter sorting the spectra according to (I-K). Linear flux den- 
sity units are used (Fa). The spectra have been normal- 
ized to a common integrated flux, and each is offset from 
the lower one by one unit. Straight segments interpolate 
through the regions where the telluric absorption was too 
strong for a satisfactory correction. 



determined structure of the latter. We w ill att empt to fur- 
ther interpret these differences in Sect. 3.1.5 . Differences 
in the width and depth of the wings of the broad H2O 
absorption bands contribute significantly to the spread 
around the mean near-IR broad band colours. CO absorp- 
tion in the H band plays a smaller role. Real deviations 
in the pseudo-continuum between 1 and 1.3 /xm occur (see 
Sect. 3.5 of LW2000) and also contribute to the spread. 
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Fig. 3. Common near-IR indices as a function of (I-K), for 
O-rich LP Vs. Lines show the properties of averaged spec- 
tra, with the three bin definitions given in the text. Plots 
(d) and (e) are flux density ratios, measured through 
the bandpasses of the HST/NICMOS/Camera2 filters 
F204M, F222M, F237M, as labeled, and converted to mag- 
nitudes (0 for Vega). The CO index in plot (f) is colour- 
corrected (see text). 



sitive the spectral signatures of molecules such as H 2 0, 
that form far in the outer atmosphere (Matsuura et al. 
2001, models of Hofmann et al. 



1998 



priv. communica- 
tion by M. Scholz), are to the instantaneous pulsation- 



Fig. 4. The seven O-rich LPVs of bin number 5 (in the 
order of Tabic [I] from the bottom upwards) . An illustration 
of the variety of near-IR properties seen at a given (I-K). 

Along the sequence of averaged spectra, the water 
vapour absorption bands deepen and widen with decreas- 
ing temperature. The measurements shown are based on 
the HST/NICMOS filters F204M and F222M: 

H 2 = -2.51og[F A (F204M)/F A (F222M)] - H 2 0(Vega) 

The first is centered at 2.04 /zm, deep in one of the H2O 
bands, the second at 2.22 /xm, at a quasi-continuum point. 
The lower envelope of the data points in Fig.^(d) illus- 
trates the mere effect of the change in colour-temperature 
for stars with shallow or inexistent H2O bands. At a given 
(I-K), star-to-star differences of up to 0.7 magnitudes are 
seen in the (204-222) index. Figure [| shows that sim- 
ilar spreads are expected in other band wings of this 
molecule, thus affecting in particular the broad band pho- 
tometry For instance, the star with the largest (H-K) (at 
(I-K)~3.8) also has extreme water bands (second spec- 
trum from the top of Fig.0). 

The CO index in Fig. g (e) measures the ratio of the 
flux densities at 2.37 /im and 2.22 /zm, based on the HST/- 
NICMOS filters F222M and F237M: 

CO = -2.51og[F A (F237M)/F A (F222M)] - CO(Vega) 

The correlation is poor, because stronger CO band 
strengths occur in redder stars and these continuum 
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changes nearly compensate for the increasing molecular 
absorption. The index shown in Fig.[|(f) includes a colour 
correction: in the above formula, F>(F237M) is replaced 
by the linear extrapolation to 2.37 /im of the flux densi- 
ties measured in the filters F165M (1.65 /im) and F222M 
(2.22 /im). With this redefinition, the flux in the CO band 
is compared to an estimate of the continuum flux at the 
same wavelength, rather than to a measurement in an off- 
set continuum filter. A strong correlation with (I-K) be- 
comes apparent. In fact, the redefined CO index combines 
the information contained in the standard CO index with 
the information contained in the colour (H-K), providing 
a better correlation with (I-K) than any of these two in- 
dices separately (the new linear correlation coefficient is 
0.82; it was 0.57 for the standard CO index and 0.77 for 
H-K). 

For comparison, we note that the LW2000 sample con- 
tains static giants with (I-K) in the range 2-3.1 and super- 
giants in the range 2.7-4.5. For all these objects, the H2O 
index of Fig.||(d) remains below 0.29. The CO indices of 
giants lie along the mean LPV lines in Figs.|](e) and (f), 
while those of supergiants lie dispersed along the upper 
envelope of the LPV distributions. 

In view of the large spread of near-IR properties found 
among O-rich LPVs, it came as a surpri se tha t the se- 
quence of average spectra described in Sect. |3.1.l| displayed 
such a regular behaviour. Regularity is essential for pop- 
ulation synthesis applications. The effects of pulsation on 
the upper atmospheres of LPVs are responsible for a large 
scatter in the properties of individual spectra, but do not 
completely wash out the underlying global effect of tem- 
perature. The sample of LW2000, for the first time, pro- 
vides enough data to recover these trends. 

3.1.3. Secondary parameters for O-rich LPV spectra 

As temperature decreases along the TP- AGB, luminosi- 
ties rise, radii increase, and masses decrease due to mass 
loss. As a consequence, stellar pulsation properties change. 
Because of the effects of initial mass, of thermal pulses and 
of LPV pulsation, stars with different pulsation properties 
can be found at a given TP- AGB temperature. Essential 
near-IR features such as the H 2 bands depend on pulsa- 
tion, even though they also correlate with effective tem- 
perature. In evolving stellar populations, the average lu- 
minosity of TP- AGB stars decreases with age and a de- 
crease of the average period is thus expected from empir- 
ical period- luminosity correlations (e.g. Feast et al. 1989, 
Wood pM Wood et al. |1999p. This example illustrates 
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that direct effects of pulsation properties on spectra at a 
given temperature will have consequences on population 
synthesis predictions. In principle, population synthesis 
models should include a gradual evolution of the pulsa- 
tion properties instead of simply considering a dichotomy 
between variable and non-variable TP- AGB stars. 

We have investigated whether empirical correlations in 
the LW2000 data between spectrophotmctric properties 
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Fig. 5. Colours and near-IR molecular indices versus pul- 
sation period, for O-rich LPVs. 



and pulsation characteristics justify the use of pulsation 
criteria rather than temperature in the definition of av- 
eraging bins. Correlations would have to be significantly 
better than those of Fig.|^, to compensate for the fact 
the theoretical predictions for the evolution of period and 
amplitude along the TP- AGB (or their empirical determi- 
nation) are even more uncertain than temperature predic- 
tions and are not included in current population synthesis 
codes. As the conclusion is negative, we only summarize 
the main empirical results here. 

In Fig.||, selected indices relevant to stellar popula- 
tion studies are plotted against period. Correlations are 
clearly present in all panels. The molecular indices show a 
larger dispersion at a given period than they do within the 
colour-temperature bins of Fig.||. The spread in (J-H) and 
(H-K) versus period is comparable to the spread seen in 
Fig. U We note that no systematic correlations are found 
between the residuals in Fig.FJ and period: for a given 
spectrophotometric index, the value and sign of the corre- 
lation coefficient changes randomly from one temperature 
bin to the next. This probably results from poor sampling 
of intrinsically dispersed relations. 

Figure^] is the equivalent of Fig.|5| for amplitude. 
Correlations between amplitude and spectrophotometric 
indices in the sample have a low level of significance, ex- 
cept for H2O (the correlation coefficient is 0.3, which for 
63 data points is significant at the 98 % level; other corre- 
lation coefficients are below 0.15). Both strong and weak 
molecular features are found at all amplitudes. We recall 
however that the smallest amplitude in the sample is of 
0.8 magnitudes and that no strong bands of water vapour 
or VO are seen in the spectra of the static or quasi-static 
giant star spectra of LW2000. Again, no systematic corre- 
lations were found between the residuals of the relations 
of Fig.|| and pulsation amplitude. 
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Fig. 6. Colours and near-IR molecular indices versus op- 
tical pulsation amplitude, for O-rich LPVs. 

With the available sample, the selection of tempera- 
ture as a primary sorting criterion for the spectra remains 
the most sensible choice. More stellar data are necessary 
to correctly define the relations between pulsation prop- 
erties and spectra within individual temperature bins. 

3.2. Carbon-rich LPVs 

Depending on initial mass and metallicity, some TP- AGB 
stars will become carbon stars through convective dredge- 
up of freshly synthesized carbon from the core, induced by 
thermal pulses (e.g. Iben & Renzini 1983| ). Carbon stars 
are essential contributors to the near-IR emission from 
AGB stars, in particular at subsolar metallicities (Persson 
et al. 1983: Frogel et al. 1990). The spectral library con- 



tains 21 complete spectra of carbon-rich LPVs. On the 
basis of theoretical spectra, Loidl et al. (2001, hereafter 



LLJ2001) showed that (R-J) and (R-H) are among the 
best effective temperature indicators for these stars. We 
will (arbitrarily) use (R-H) here. The far-IR source RLep 
is so red that it must be treated separately. RLep ex- 
cluded, the range of (R-H) values in the sample does not 
justify more than two or three averaging bins. 

The suggested temperature sequence for C-rich stars 
contains 3 averages of 6 spectra each (see Tablc^), fol- 
lowed by the individual spectrum of R Lep near maximum 
light, and finally the average of the 3 RLep spectra near 
minimum light (the optical spectrum in the 4th bin is an 
average of bins 3 and 5, because no data was obtained). 
It is presented in Fig. f[| The averaging procedure is as for 
O-rich LPVs. The spectrophotometric properties of the 
data bins are shown in Fig.||. The C2 index measures the 
strength of the bandhead at 1.77/im. To avoid the spec- 
tral region of strong telluric absorption corrections and 
of contamination by CN, the "continuum" flux density 
is measured over the narrow range 1.752-1.762 /im and 



the b and flux density over 1.768-1.782 //m (Alvarez et al. 
|2000D ; the index is given in magnitudes. The models of 
Gautschy ( 2001 ) demonstrate that the C2 bandhead is 
sensitive to the C/O ratio, the 12 C/ 13 C ratio and micro- 
turbulence. This contributes to the dispersion in panel (e). 
More data is needed to assess the significance of the trend 
in panels (e) and (f). 




#^V-k^ fw. 



CN j.jM 



/ --|WvV»V»'V' ' '"'~' r 

1K1 



0.5 1 1.5 2 2.5 

Wavelength (/im) 

Fig. 7. A temperature sequence of carbon star spectra. 
Linear flux density units are used (Fa , after normalization 
to a common total flux, with offsets of 1.1). In the regions 
irremediably affected by telluric absorption, around 1.35 
and 1.9 /im, we have drawn segments by hand that approx- 
imately follow the theoretical shape of molecular bands 
(LLJ2001). The carrier of the molecular band around 
1.55 /im is uncertain (Gautschy 2001 suggests a combi- 
nation of C 2 H 2 and HCN). 



Alternative criteria in the construction of averaging 
bins for carbon stars are C/O or 12 C/ 13 C, because their 
effects on the spectral signatures are significant. Both vary 
along the TP-AGB as a consequence of third dredge-up 
episodes and envelope burning, in a way that is becoming 
predictable (Renzini & Voli 1981, Forestini & Charbonncl 
1994 Marigo |200lD . Binning with respect to C/O can 
be done on the basis of estimates of LLJ2001. It is best 
to isolate the S/C star BH Cru in one bin, to combine 
T Cae and RU Pup data, and to group Y Hya and S Cen. 
Unfortunately, current atmosphere models do not allow us 
to estimate the ratio for stars as cool as RLep. The em- 
pirical sequence of increasing C/O ratios is characterised 
by increasing C2 band depths (at 1.77 /im, around 1.2/im 
and below 5700 A), and by decreasing apparent strengths 
of the CO bandheads around 1.6 /im and 2.3/im as a re- 
sult of heavier contamination by CN and C2 (cf. Fig. 4 
of LLJ2001). The C/O and 12 C/ 13 C sequences based on 
the LLJ2001 estimates are similar (though in reverse or- 
der). We note however that a degeneracy with the effects 
of turbulent velocity makes 12 C/ 13 C estimates uncertain. 
The inverse correlation between C/O and 12 C/ 13 C found 
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Fig. 8. Spectrophotometric properties of carbon stars as 
a function of the bin number based on (R-H) . The indices 
for instantaneous spectra are shown as crosses; the solid 
lines connect those of the averaged spectra (bin 4 contains 
only 1 spectrum). 

Table 2. Contents of the C-rich temperature bins 



Bin 


Name 


Date 


C/O a 


Period 


8V 


(R-H) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


TCae 


Dec. 95 


1.05 


156. 


1.8 


3.55 


1 


TCae 


Jan. 96 


1.05 


156. 


1.8 


3.58 


1 


SCen 


Jan. 96 


1.40 


65. 


1.5 


3.78 


1 


TCae 


Mar. 96 


1.05 


156. 


1.8 


3.83 


1 


BHCru 


Jan. 96 


1.01 


421 


2.8 


3.84 


1 


YHya 


Dec. 95 


1.40 


302.8 


3.7 


3.89 


2 


RUPup 


Dec. 95 


1.10 


425. 


1.9 


3.98 


2 


RUPup 


Jun. 95 


1.10 


425. 


1.9 


4.00 


2 


YHya 


May 96 


1.40 


302.8 


3.7 


4.04 


2 


YHya 


Jan. 96 


1.40 


302.8 


3.7 


4.07 


2 


YHya 


Mar. 96 


1.40 


302.8 


3.7 


4.10 


2 


YHya 


Jul. 96 


1.40 


302.8 


3.7 


4.20 


3 


RUPup 


Jan. 96 


1.10 


425. 


1.9 


4.22 


3 


RUPup 


Mar. 96 


1.10 


425. 


1.9 


4.25 


3 


YHya 


Jun. 95 


1.40 


302.8 


3.7 


4.26 


3 


BHCru 


May 96 


1.01 


421. 


2.8 


4.40 


3 


RUPup 


May 96 


1.10 


425. 


1.9 


4.47 


3 


BHCru 


May 96 


1.01 


421. 


2.8 


4.78 


4 


RLep 


Nov. 97 




427.1 


6.2 


5.08 


5 


RLep 


Dec. 95 




427.1 


6.2 


5.58 


5 


RLep 


Jan. 96 




427.1 


6.2 


5.94 


5 


RLep 


Mar. 96 




427.1 


6.2 


6.16 



Note to TaMe|2j: 

a: Atmospheric number ratio of carbon to oxygen, when avail- 
able, according to the estimates of LLJ2001. 



in the LW2000 sample, and values of 12 C/ 13 C as low as 
3.5 are not typical from the stellar evolution point of view. 

Pulsation has small effects on the spectrum below 
2.5 fim, except for the coolest carbon stars, where it allows 



a dee p band to form around 1.55 /xm (Fig.^ and Gautschy 

EootD . 



3.3. Dust-enshrouded LPVs 

OH/IR stars and their carbon rich counterparts are pro- 
duced on the upper AGB when the dust produced in the 
heavy winds of an LPV builds a thick enough envelope to 
hide the central star at optical wavelengths. At a given 
bolometric luminosity, their contribution to the near-IR 
light is smaller than for unobscured stars but their mid- 
IR emission is significant, an effect that may be useful in 
searc hes for intermediate age stellar populations (Bressan 
et al. 



1998|) 



More than the effective temperature of the 
central star, the amount of dust (which is related to the 
mass loss and, indirectly, also to the effective temperature) 
determines the near-IR energy output. 

To represent the optical and near-IR emission of dust- 
enshrouded sources, we suggest to use reddened versions 
of the average spectra provided in this paper. Although 
circumstellar dust might be of a different nature than 
the average interstellar dust grains, the LW2000 data sug- 
ges t tha t the Milky Way extinction law (e.g. Cardelli et 

3.1) can be used to first or- 



al. |1989| , with Av/E(B-V)= 
der for that purpose. The library of LW2000 contains a 
handful of spectra of OH/IR sources. In Fig.|, AFGL1686 
and WXPsc are compared with reddened versions of the 
coolest average spectrum of Sect. 3.1.1, which has similar 
bands of H2O, TiO and VO. The agreement is rather good 
(and even better agreements can be found with reddened 
versions of selected individual spectra of the initial spec- 
tral library). The Milky Way extinction law creates no 
obvious discrepancy. More data will be needed to inves- 
tigate this question further, and to perform similar tests 
with carbon rich stars. 

4. Using the spectra in population synthesis codes 

For practical use in connection with theoretical stellar 
evolution tracks, the empirical spectra must be properly 
normalized and a quantitative effective temperature scale 
must be adopted. These aspects are discussed in Sect. 4.1 
and Sect.[h^ below. Of course, the stellar evolution tracks 
must also provide the transitions between O-rich and C- 
rich atmospheres, and between optically visible sources 
and dust-enshrouded infrared sources. The discussion of 
these aspects of theoretical stellar evolution lies outside 
the scope of this paper. 

4.1. Bolometric corrections 

Bolometric corrections describe the relation between the 
spectral flux densities and the bolometric luminosity Lboi- 
As the LW2000 data encompass most of energy emitted by 
the LPVs, only small corrections are necessary to convert 
the measured integrated flux into a bolometric flux. For 
oxygen rich stars, Alvarez et al. ( [2000 ) computed these 
corrections by extrapolating the LW2000 data with the 
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Fig. 9. (a): Average of the April 1995 and May 1996 
spectra of AFGL 1686 = IRAS 14068-0730 (dotted ) and 



reddened version of the spectrum of bin 9 of Sect. |3. 1.1 
using the Milky Way extinction law with Av=3 (solid), 
(b): Spectrum of WXPsc = IRAS 01037+1219 (dotted) 
and reddened spectrum of bin 9 with Av=18 (solid). Fa 
units are arbitrary. 
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best-fitting cool M giant models. They found that on av- 
erage about 20 % of the light is emitted outside the range 
of the available data. In practice, because blackbodies are 
poor representations of the infrared emission of cool stars 
(Fluk s et al. 1994 , Aringer et al. 1996 , J0rgensen et al. 
2001 ), we suggest to simply extrapolate the data linearly 
to a nil flux at 3500 A on the blue side and to 50 000 A 
on the red side. These extrapolation limits are chosen 
to reproduce the results of Alvarez et al., as shown in 
Fig.[l0|(a). The long wavelength extrapolation also pro- 
vides agreement with the integrated fluxes measured long- 
ward of 2.5^m for the six ISO satellite spectra of oxygen 
rich LPVs we have examined^], to within random devia- 
tions smaller than 10% (i.e. a few percent on Lboi)- 

The infrared emission of carbon stars also differ s sig- 
nifica ntly f rom a blackbody (e.g. Goebel et al. 197£ , Aoki 
et al. 1996 ). According to a subset of the models described 
in LLJ2001, and found to provide satisfactory agreement 
with ISO data (J0rgensen et al. |2000|) , the wavelength 
range of the LW2000 data again comprises 80 % of the 
flux, and only a few percent of the missed flux comes 
from A < 5000 A. The total flux of ISO spectra]] longward 
of 2.5 fim is equivalent to the suggested linear extrapo- 
lation, with object-to-object differences that amount to 
variations of up to 10 % in the total fluxes. The absence of 
ISO counterparts for the LW2000 stars makes it impossi- 
ble to test directly whether the extrapolation is still valid 



b SVPeg, TSge, V1351Cyg, V584Apl, VYCas 
and CEAnd; available from the public domain page 
of Kerscfibaum, Loidl, Hron, Kerber, and Rauch: 
www.astro.univie.ac. at/ ~fzi/AGB/agb2pn. html. 

7 TXPsc, V460 Cyg and SCep; retrieved from the Vienna 
web site mentioned above. 



Fig. 10. Bolometric corrections, obtained from the 
LW2000 data and the mean spectra, with the extrapo- 
lation method described in the text. Crosses and open 
squares: individual stellar spectra for O-rich and C- 
rich LPVs; solid symbols: sequences of averaged spec- 
tra. Dotted line in (a): values obtained by Alvarez et al. 

( pooo| ). 



for the coolest carbon star bin, but the near-IR similarity 
between RLep and the ISO target SCep (see Goebel et 
al. 1978 and LW2000) suggests that the linear extrapola- 
tion method doesn't underestimate the Lboi by more than 
30 % even in that case. 

Figure [To] shows that the bolometric corrections 
Mboi— I and Mboi— J for carbon stars lie below those of 
oxygen rich stars. The same would be true for Mboi- H. 
The explanation lies in the location of the main molecu- 
lar bands. In O-rich stars, water bands tend to lie at the 
edge of the standard photometric windows, while in car- 
bon stars the C2 and CN bands reduce the flux density 
right through these windows. In comparison, molecular 
absorption is relatively small in the K window for both 
chemical classes. 

The bolometric corrections above do not apply to AGB 
stars in opaq ue d ust shells. Applying reddening as sug- 
gested in Sect. 3.3 automatically accounts for the reduced 
optical and near-IR fluxes of these objects. The determi- 
nation of bolometric corrections is thus directly dependent 
on the rapidly evolving circumstellar optical depth. 
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4.2. Temperature scales 



Table 3. Properties of the averaged spectra 



The second step required in the assignment of a spectrum 
to a particular stage along a stellar evolution track is the 
determination of an effective temperature (T e gp) scale for 
the available spectra. A straightforward way of doing this 
is to rely on theoretical or semi-empirical scales (e.g. an- 
gular diameter measurements) . Although we provide rela- 
tions of that kind here, we emphasize that with current 
stellar models all attempts to relate temperatures assigned 
this way with the T G ff of a particular set of stellar evolu- 
tion tracks are uncertain. More arguments to this impor- 



O-rich spectra 



tant point are given in Sect. 5.1 



Bessell et al. ( 1989 b ) provide temperature scales based 
on (I-K). The scales are metallicity-dependent. We adopt 
solar mctallicity here and refer to Sect.|| for a discus- 
sion of metallicity effects. Synthetic spectra for static and 
pulsating AGB stars have been kindly made available to 
us by M. Scholz. The static models and t heir colours are 
described in Bessell et al. ( 1989 b| , 1991), the pulsating 



ones in Bessell et al. (1996; their E, Z and D series) and 
Hofmann et al. ( 1998| ; their L, M, P and Q series). A de- 
tailed comparison between the models and the data lies 
far beyond the scope of this paper. Although many ob- 
served tendencies are reproduced by the models at least 
in a statistical sense, the fits to individual spectra are not 
perfect, and from our preliminary analysis we estimate 
that they do not provide a more reliable absolute effective 
temperature than, for instance, a colour such as (I-K). 




Fig. 11. Effective temperature scales based on synthetic 
spectra. The lines follow static model series. Dashed 
(Bessell et al. 1989b): 1M Q sequences with, from top 
to bottom (10 4 L Q ; Z Q ), ( 210 3 L Q ; Z Q ), (1O 4 L ; Z /2). 
Solid (Bessell et al. 1991): the bundle of 4 similar lines 
follow the AGB locus at Z Q for 1, 1.5, 2.5 and 5M , and 
the bluer line is for 1 M© at Z /2. Symbols (Bessell et 
al. 1996, Hofmann et al. 1998| ): dynamical model stars at 



Bin 


I-K 


J-K (J-Ks) a 


BC(I) 


T 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


1.96 


0.96 (1.01) 


0.65 


3930 K 


2 


2.40 


1.08 (1.13) 


0.43 


3585 K 


3 


2.76 


1.17 (1.22) 


0.23 


3440 K 


4 


3.10 


1.18 (1.23) 


-0.04 


3355 K 


5 


3.45 


1.29 (1.34) 


-0.28 


3285 K 


6 


3.79 


1.26 (1.31) 


-0.57 


3235 K 


7 


4.28 


1 on / 1 o /i \ 
1.20 (1.24) 


-1.03 


3175 K 


8 


4.88 


1.38 (1.43) 


-1.45 


3055 K 


9 


5.60 


1.50 (1.55) 


-2.01 


2340 K 




C- 


rich spectra 






Bin 


R-H 


J-K (J-Ks) a 


BC(I) 


T 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


3.74 


1.43 (1.52) 


0.02 


3200 K 


2 


4.06 


1.63 (1.71) 


-0.06 


3000 K 


3 


4.38 


1.65 (1.74) 


-0.30 


2800 K 


4 


5.08 


2.15 (2.23) 


-0.97: 


2400 K 


5 


5.86 


3.03 (3.08) 


-1.96: 


2000 K 



Notes to TaWeg: 

a : (J-Ks ) is given for the DENIS filter passbands (Epchtein et 



al. 1997), kindly provided by M. Schultheis. 
A colon indicates uncertain values, where the far-IR emission 
is likely to be underestimated by the adopted extrapolation 
procedure. 



For the interpretation of Col. 5, see Sect. 4.2 and 5.1 



Figure [Tl| shows the relation between (I-K) and the ef- 
fective temperature for solar metallicity models and for 
two series of static models at subsolar metallicity. The 
effective temperature is a theoretical quantity, defined us- 
ing Stefan's law and the radius at which the Rosseland 
mean optical depth of the star reaches unity. While each 
of the dashed lines is a sequence of models at con- 
stant mass (1M ) and luminosity, the solid lines follow 
the luminosity-temperature relations of AGB evolution- 
ary tracks as parametrized by Bessell et al. (1991). The 
symbols shown for the pulsating models are instantaneous 
values, corresponding to a particular cycle and pulsation 
phase. It is seen that the dynamical models follow a se- 
quence that agrees with the IMq, high luminosity (i.e. 
low gravity) static solar metallicity sequence of Bessell et 
al. ( 1989 b| ). This locus provides the temperatures of the 
averaged empirical spectra listed in Col. 5 of Table||. A 
lower effective temperature would be assigned to a given 
spectrum if the subsolar theoretical scales or the solar ones 
of Bessell et al. (1991) were used instead. The compari- 



phase 0.5 (triangles) or 0.8 or 0. (squares) 



son of the range in (I-K) in Table ^ and along the models 
in Fig. [ll] shows that the sequence of average LPV spec- 
tra adequately covers the temperature range required for 
AGB synthesis. 

For carbon stars with effective temperatures above 
~ 2600 K, we adopt the relation between (R-H) and T e fi 
derived by LLJ2001. The values for the binned spectra of 
Fig. ^ are given in Table ||. The two coolest bins of the 
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sequence contain spectra of RLep, a Mira-type variable 
with a period of 427 days. A T off of 2058 ± 180 K has been 
obtained for this star from interferometric angul ar siz e 
measurements in October, 1995 (van Belle et al., 1997 ). 



The data in our temperature bin 5 were taken between 
December, 1995 and early March, 1996. The 3 individual 
spectra are very similar; the comparison of their energy 
distributions suggests that the minimum temperature was 
reached in early January, 1996, i.e. only 0.2 cycles after 
the interferometric measurement. We therefore assign the 
spectrum of bin 5 an effective temperature of ~ 2000 K. 
The (V-K )-T e ff relation for Carbon Miras of van Belle et 
al. (1997) predicts (V-K) ~ 10 at this temperature, which 
is in good agreement with our measurement (after extrap- 
olation of the observed spectrum to the short wavelength 
cut-off of the standard V band filters). Based on this se- 
quence, T e fj ~ 2400 K is appropriate for bin 4. The range 
of (J-K) values in the average spectra is in excellent agree- 
ment with observations in complete carbon star samples 
(Frogel & Elias |1988| , Glass et al 
Weinberg & Nikolaev |200l|) . 



1995, Cioni et al. 200C 



5. Discussion 

5.1. Temperature bins and temperature calibrations 

In the above sections, we have shown that wide baseline 
colour temperatures provide the most sensible and prac- 
tical classification of the LPV spectra. When using these 
binned spectra to synthesize cluster or galaxy spectra, one 
assumes that the energy- weighted mean spectrum of an in- 
dividual LPV ( averaged over the < 10 3 days of the LPV 
pulsation cycle) is similar to the average of many individ- 
ual stellar spectra obtained from stars of various masses, 
amplitudes and phases, but with a common colour temper- 
ature. As we found no significant correlations, in any of our 
temperature bins, between the deviations from the mean 
spectral properties and parameters such as amplitude or 
period, this is a reasonable choice. More observations will 
be needed before other options can be considered. 

With the above assumption, the effective tempera- 
ture along a TP- AGB evolution track alone determines 
which spectrum should be used to represent the corre- 
sponding emission (it is assumed that the track also indi- 
cates when the use of LPV spectra instead of static giants 
becomes necessary, and when a star becomes C-rich or 
dust-enshrouded). In order to avoid misinterpretations of 
the temperature scales given in Table^, it is important 
to recall how evolutionary TP- AGB tracks and individual 
TP- AGB stars behave in the HR diagram. 

Even a perfectly coeval stellar population produces a 
broad asymptotic giant branch in the HR diagram and 
in colour-magnitude diagrams, as a result of the thermal 
pulses (helium shell flashes) that occur on timescales of 
~ 10 4 yrs and of the long period pulsational variability 
on timescales of 10 2 — 10 3 days. Figure [l^ shows theoret- 
ical illustrations of these effects, and Fig.|l^ provides an 
empirical counterpart. 




Fig. 12. Left : Synthetic TP- AGB track for the evolu- 
tion of a star with an initial mass of 2 M Q and a final 
mass of 0.6 Mq. Right : 1 Mq pulsating star models from 
Hofmann et al. (1998). The P, M and O series are shown 
respectively as squares, circles and triangles. Properties 
at maximum and minimum light in several successive cy- 
cles are shown (open symbols) , together with those of the 
static parent star (solid). 



The TP- AGB track on Fig.|| (left) was computed 
using so -calle d "synthetic evolution" t echni ques (Iben & 
Rcnzini 1983, Groenewegen & de Jong 1993, Wagenhuber 
Groenewegen 1998, Marigo et al. 



1998, Mouhcine & 



Langon [2002 ). The values of the various model parame- 
ters for this particular track are not relevant here. Tiny 
offsets in many stellar parameters will ensure that thermal 
pulses do not occur in phase for all stars. Most stars will 
be found around the quiescent luminosity, and 10 — 30% 
at the significantly lower luminosities that follow thermal 
pulses (Boothroyd & Sackmann 1988, Groenewegen & de 
Jong 1993, Marigo et al. 1999| ). The locus of any set of 



TP-AGB tracks is based on its own particular definition 
of the T e ff (see Scholz & Takeda |1987 for a discussion), 
through rather simple outer boundary conditions for such 
cool and extended stellar atmospheres. In addition, as al- 
ready mentioned in Sect. [l|, these temperatures are those 
of static model stars. 

The right hand plot of Fig. [l^ shows the predicted ef- 
fect of Mira-type pulsation for the three most recent mod- 



els described by Hofmann et al. (1998). All models pulsate 
with periods of the order of 330 days and assume a mass 
of IMq. The two models at relatively lower luminosity 
(series P and M of the authors) are fundamental mode 
pulsators, while the third one (O series) pulsates in the 
first overtone. All models behave differently in the HR di- 
agram. Model M shows that pulsation is able to shift a 
star systematically to larger radii and lower effective tem- 
peratures, when compared to its static parent. This trend 
is also seen when direct angular diameter measurements 
of pulsating and static late-type giants are compared (van 
Belle et al. 1997). However, the amplitude of this effect at 
a given point of an evolutionary track cannot be predicted 
reliably because neither the evolution of LPV pulsation 
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Fig. 13. Tentative location of the LW2000 spectra on the 
theoretical HR diagram (for stars with known periods). 
T eB is derived from (I-K), log(L) from the period and the 
instantaneous bolometric correction. Symbols distinguish 
individual stars; multiple observations at various phases 
are connected. Smoot hed A GB tracks at Z Q are overlaid 
(Mouhcine & Langon fc002[ at 0.8 M , only the RGB is 
shown) . 

modes, nor pulsation amplitudes, nor the actual effect of 
these quantities on the spectra are understood. 

In Fig. p], the O-rich stars of the LW2000 library have 
been located tentatively in the theoretical HR diagram, in 
order to estimate the empirical width of an upper AGB 
sequence. Each star with a known period has been as- 
signed a K band luminosity (Lk) using the empirically 
well-defined period - Lk relation for Miras by Hughes & 
Wood ( 1990| ). Note that those stars which pulsate in over- 
tone modes should in fact be assigned higher luminosi- 
ties (see LW2000 for potential candidates). Lk variations 
have been neglected in the figure. The changes in the bolo- 
metric luminosity represent the changes in the instanta- 
neous b olom etric corrections measured on the spectra as 
in Sect. 



4.1 



The assigned temperatures are derived from 
(I-K) as in Sect. 4.2. It is seen that pulsation can widen the 
AGB to AT = 500- 1000 K or A(I - K) = 1-2 mag. This 
range is consistent with the predicted ranges shown in 
Fig.|l2|. If pulsation was neglected, such a range of colours 
could be mistaken for a range of several M Q in initial stel- 
lar masses. Alternatively, it could be attributed to a spread 
in metallicity. TP- AGB tracks are typically shifted to hot- 
ter T cff by 0.3 dex when going from Z=0.02 to Z=0.008. 

The tracks overlaid on Fig.[l3] are smoothed (many 
population synthesis codes require that thermal pulses be 
smoothed out, to make interpolations between tracks or 



the construction of isochrones possible). The smoothed lu- 
minosities are the result of energy conservation: smooth 
tracks run at luminosities slightly lower than the qui- 
escent, pre-flash luminosities, in order to account for 
the post-flash luminosity dips. More options exist in the 
smoothing schemes for T c g , as the ideal weighting should 
use the energy output in the particular wavelength range 
one intends to make predictions for. The differences be- 
tween various schemes will however remain small com- 
pared to the systematic effects of T e ff definitions or LPV 
pulsation mentioned above. 

In summary, theoretical effective temperatures are not 
directly comparable to those estimated for individual em- 
pirical spectra of LP Vs. Acquiring enough data to obtain 
representative average spectra is only one step, in which 
the LW2000 sample has allowed significant progress. But 
until consistent models for the structure and spectrum of 
pulsating AGB stars exist, the link between the theoretical 
T e ff scale and any semi-empirical one remains uncertain. It 
must be calibrated a posteriori. Ideal observational targets 
for the temperature calibrations are massive stellar clus- 
ters with known metallicities and ages (from optical spec- 
troscopy), such as found in merger galaxies for instance 
(Maraston et al. |200l| , Mouhcine et al. |2002|) . The rele- 
vant ages lie between 10 8 and 2 10 9 yrs. In individual clus- 
ters containing less than ~ 10 5 stars, the stochastic effects 
due to the intrinsically small number of bright AGB stars 
affect the comparison significantly (Langon & Mouhcine 
|2000D . 



5.2. Period and amplitude distributions 

As recalled in Sect.||, the stars of the LW2000 sample 
were selected to provide a relatively uniform coverage of 
the period-amplitude plane. Those used to construct the 
averages of this paper cover periods of 100 < P < 450 and 
optical amplitudes of 0.8 < <5V < 7 magnitudes. 

Systematic photometric observations of large samples 
of stars have shown that the majority of the TP- AGB stars 
are variable (Wood et al. 1999, Alcock et al. 2000, Cioni 
et al. [2001 ). However, period and amplitude distributions 
are expected to depend on the age of the observed pop- 
ulation (Vassiliadis & Wood 1993). Large amplitude pul- 
sation characterizes only a subpopulation of the TP- AGB 
stars. Here, we discuss whether or not it is appropriate 
to use the average LPV spectra of this paper for all the 
TP- AGB stars of a synthetic population. 

Period distributions in the Large Magellanic Cloud 
(LMC) or in the Galaxy extend from less than 20 days to 
more than 10 3 days. The range of periods of our LPV sam- 
ple corresponds to those of optically visible Miras and of 
relatively long period, high luminosity semi-regular stars, 
as seen in Fig. 1 of Wood ( |2000[ ). Both in the LMC and 
in the Galactic Bulge, semi-regular variables with peri- 
ods shorter than 100 days outnumber Miras (Wood et al. 
1999| , Alard et al. ^00 1[ Alcock et al. |2000| ). The mean pe- 
riod and the mean amplitude of our sample are thus larger 
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than the mean period of optically visible TP- AGB stars in 
populations that have been forming stars for a long time 0.4 
(i.e. since the epoch of galaxy formation). 



What is relevant to the synthesis of integrated galaxy 
or cluster spectra is how the period and amplitude dis- 
tribution of the available stellar sample compares with 
the distribution for those particular TP-AGB stars that 
contribute most strongly to the light (in the near-IR, i.e. 
the only spectral range of the LW2000 data where the 
TP-AGB contribution is significant). Population synthesis 
models have shown that these are stars with initial masses 
between about 1.8 and 3M Q (Mouhcine & Langon 2002; 
see also Renzini 1992j , Girardi & Bertelli 1998, Maraston 
1998] ). They reach the TP-AGB with ages of about 0.1 to 
1.3 Gyr and spend of the order of 10 6 yrs in that phase (ex- 
act ranges are model and metallicity dependent). In older 
populations, the stars of the red giant branch are so nu- 
merous that they provide most of the near-IR light. Which 
spectra are used for the remaining insignificant TP-AGB 
stars stops being critical then. 

What do we know about the pulsational properties 
of stars with initial masses between 1.8 and 3M Q ? A 
global picture has emerged from joint studies of empir- 
ical period-luminosity distributions and theoretical pul- 
sation models (Fox fc Wood 1982, Xiong et al 



1998 



Vassiliadis k Wood |l993| , Wagenhuber & Tuchman |1996|) 



In the period-luminosity diagram, a single evolutionary 
track runs "stepwise" accross pulsation mode sequences, 
from short to long periods, with a slowly increasing lu- 
minosity (apart from the temporary changes due to the 
thermal pulses): a star becomes unstable to one of the 
overtone pulsation modes, follows that period-luminosity 
sequence until it switches to a lower overtone, and so on. 
The large mass loss associated with Mira-type pulsation 
finally makes the star loose its envelope. Stars with more 
massive initial masses arrive in the unstable strips of the 
period-luminosity diagram with larger luminosities. This 
is supported by the observation that Milky Way Miras 
with very long periods are associated with the thin disk, 
while shorter period Miras are seen in the thick disk and 



in globular clusters (Feast |1963 Habing 199£). The stars 
with masses between 1.8 and 3M Q are thus expected to 
pulsate with larger mean periods than a population in 
which most stars are old. 

Figure [l4] shows predicted period distributions for a 
stars with an initial mass of 2 Mq, based on the TP-AGB 



tracks of Mouhcine & Langon ( [2002 ), c ombin ed with the 
period-mass-radius relations of Wood (|1990D for funda- 
mental mode and of Barthes & Luri ( 2001 ) for the first 
overtone. Systematic effects of the initial mass are small 
between 1.8 and 3M Q . The plotted distribution is rep- 
resentative of the initial masses most relevant to observ- 
able integrated properties of stellar populations. The up- 
per graph is a number distribution, the lower gives the 
relative contributions of stars with various periods to the 
luminosity. The relative number of fundamental mode and 
overtone pulsators is very uncertain. In the LMC period- 
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Fig. 14. Predicted period distributions for stars with ini- 
tial masses around 2M Q (Z=Z©). Fundamental mode pul- 
sation is assumed for the solid line, first overtone for the 
dashed line. Each distribution is normalized individually 
to a sum of one. 



luminosity diagram (Wood et al. 1999), about 1/3 of the 
LPVs in the relevant, high luminosity range are assigned 
fundamental mode (unfortunately the LMC scarcely pop- 
ulates the upper part of the diagram, due to the com- 
bined effects of the stellar initial mass function and the 
star formation history). Fig.|l4| shows that the range of 
periods available to us for the construction of average 
spectra is indeed very appropriate for the fundamental 
mode pulsators. A similar conclusion is reached when 
other tracks in the literature are used (Vassiliadi s fc W ood 
1993j Wagenhuber & Tuchman |l99^ , Gautschy |l999| ) . 

The periods of our spectroscopic sample are however 
larger than those of the overtone pulsators. According to 
the data of Wood et al. ( 1999] ), the overtone pulsators 
are essentially all semi-regular variables, with small am- 
plitudes compared to those of the Mira-like fundamental 
mode pulsators (semi-regular LPVs are found on the fun- 
damental mode sequence, but at lower luminosities and 
shorter periods than those relevant here, as seen in the 



light curves given in Wood et al. 1999). Many of the over- 
tone pulsators should be well represented with spectra of 
static giants. Thus, the relative importance of the LPVs 
in synthetic populations can be varied by modifying the 
blue edge of the instability strip. 

Current pulsation theory does not put strong con- 
straints on this edge and on the transition from overtone 
to fundamental mode pulsation. In order to evaluate the 
effects of such adjustments, we have compared synthetic 
spectra of intermediate age populations computed with 
two assumptions: (i) our average LPV spectra are appro- 
priate for all TP-AGB stars or (ii) the LPV spectra only 
apply in the Mira instability strip. For (ii), the blue edge 
of the Mira instability strip was projected from the empir- 
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ical period-luminosity plane into the HR diagr am as sum- 
ing fundamental mode pulsation models (Wood 199Cl| ), and 
static giant stars were used up to that point. The differ- 
ences were much smaller than the uncertainties currently 
associated with, for instance, the effective temperature 
scales of the spectra (changes were smaller than 3 % in 
all the quantities we measured). 

In summary, the simplifying assumption that the aver- 
age LPV spectra of this paper apply to all TP- AGB stars 
is reasonable, when the aim is the synthesis of integrated 
spectra of stellar populations. Better theoretical predic- 
tions for the onset and evolution of period and amplitude 
of TP- AGB stars with various initial masses will be needed 
to reduce the uncertainties inherent to this assumption. In 
the study of resolved populations that include AGB stars 
with low initial masses, it is likely that the spectra of static 
giant stars should be used for some as yet poorly defined 
first part of the TP- AGB. 

5.3. Hydrogen recombination lines 

The only conspicuous near-IR emission line seen in LPVs 
is Paschen /3, at 1.28 /im. Emission equivalent widths of up 
to 6 A are observed in individual spectra of the LW2000 
sample. Pa/3 is seen in most of the averaged O-rich spectra 
of Sect. H The emission equivalent width is of the order of 
1 A in 6 of the 9 temperature bins. None of the C star spec- 
tra of LW2000 shows significant Hydrogen recombination 
lines. This is mostly due to molecular line crowding, which 
is more severe in C-rich spectra around the wavelengths of 
the hydrogen lines. In the maximum light spectrum of the 
S/C type star BHCru, Pa/3 is detected with an equiva- 
lent width of the order of 1.5 A (the resolution of the data 
defines the pseudo-continuum). Paschen /3 is not seen in 



the averaged C star spectra of Sect. 3.2 



According to the evolutionary tracks and popu lation 
synthesis calculations of Mouhcine & Langon ( [2002 ) , TP- 
AGB stars contribute less than 40 % of the J band light in 
coeval populations of intermediate ages. If the Pa/3 equiv- 
alent widths of the average LPV spectra are representative 
of the TP- AGB, those of integrated spectra of stellar pop- 
ulations will not exceed fractions of an Angstrom. This is 
small compared to the effect of any significant subpopula- 
tion of young ionizing stars (Leitherer et al. 1999). 



6. A recipe for spectral synthesis at various 
metallicities 

When metallicity is varied, several effects add up to pro- 
duce differences in the predicted contributions of upper 
AGB stars to the integrated light of a stellar popula- 
tion. First, the relation between effective temperature and 
spectrophotometric properties changes. For cool stars, this 
is due to the importance of molecular opacities. An il- 
lustration can be found, e.g., in Fig. 10 of Hauschildt et 



licity, a bluer spectrum with weaker molecular bands 
must be associated with each given point of the theoret- 
ical HR diagram. Second, the AGB evolutionary tracks 
themselves are globally shifted to lower effective tempera- 
tures. This effect acts in the same direction as the first 
one: signatures of cooler AGB stars will be more pro- 
nounced at higher metallicities. The third effect is the 
metallicity-dependence of the mass loss, of the AGB life- 
times and of the efficiency of the production of carbon 
stars. Discussions of these lat ter as pects can be foun d, e.g ., 
in Groenewegen & de Jong ( 1993; ), Marigo et al. ( 1998 ), 
Mouhcine & Langon fl2002[ ). 

The metallicites of the LW2000 spectra are not known. 
As the observed stars are located in the solar neighbour- 
hood, they are expected to be representative of stars in ap- 
proximately solar metallicity environments. The range of 
ages (or initial masses) in the data is likely to correspond 
to a range of metallicities, as determined by the chemical 
evolution of our neighbourhood. A contamination of the 
sample by Population II stars is also likely (LW2000). The 
calibrations of Ramirez et al. 200C and Frogel et al. |2001 



that relate [Fe/H] to the equivalent widths of near-IR ab- 
sorption features and colours based on Milky Way cluster 
red giant branch stars, would lead to dispersed, on aver- 
age slightly subsolar metallicities. However, the calibra- 
tion samples don't include LPVs, and the observed vari- 
ability of the relevant spectral signatures in the LW2000 
sample makes the results uncertain. The LW2000 have also 
been c ompare d to the synthetic red giant spect ra of Bessell 
et al. ( |1989 aD and of Hauschildt et al. fll999| ). Both sets 
provide decent (though far from perfect) agreement with 
the observed giant star spectra and at least so me o f the 
LPV spectra. In the models of Hauschildt et al. ( 1999 ), the 
effect of moving from [M/H]=0 to -0.3 is so strong at low 
temperatures that the low metallicity model spectra are 
very different from any observed one we have seen, in the 
LW2000 data or in the literature. With these models, one 
would conclude that all the LW2000 stars have quasi-solar 
metallicities, but our impression is that the effect of Z in 
the models is too strong. In the model set of Bessell et al. 
( |l989~al ) moving from log(Z/Z Q )=0 to —1 has a relatively 
small effect on the spectra (at a given T e ff, mass and sur- 
face gravity) , indicating that the observed spectra may be 
appropriate over that range of metallicities. The compari- 
son between properties of stars in the Galactic Bulge and 
in the Field supports a modest evolution of intrinsic spec- 
tral s ignat ures over a useful range of metallicites (Frogel 
et al. |1984D . 

In view of this situation, we suggest to restrict the in- 
clusion of metallicity effects in spectral population synthe- 
sis models to those that are intrinsic to stellar evolutionary 
tracks, and to the use of a metallicity dependent T c g scale 
for the mean O-spectra. For instance, one may use the 
relations between (I-K) and T c ff given for various metal- 
licities by Bessell et al. 1991. The sensitivity of carbon 



al. ( 1999 ), where static giant model spectra for [M/H]=0 
and -0.3 are compared. When working at lower metal- 



star spectra to metallicity is small. R. Gautschy (2001) has 
computed model spectra for hydrostatic spherical carbon 
stars between log(Z/Z Q )=0 and -0.3. The models show 
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changes in Z have small effects compared to changes in 
the effective temperature, the G/O ratio or the 12 C/ 13 C 
ratio. The latter are determined by the current evolution- 
ary status of a star, at any initial mctallicity. For carbon 
stars, we suggest to use the average spectra of this paper 
with the same effective temperature scale at all Z. 



7. Summary and conclusions 

This paper provides convenient and representative data 
for stars on the upper asymptotic giant branch, to be used 
in the studies of cool stellar populations of galaxies. 

We have used the data of LW2000 to construct se- 
quences of average spectra for O-rich and C-rich LPVs. 
The O-rich sequence is based on (I-K), which is taken as 
a first order effective temperature indicator. Despite the 
large dispersion between (I-K) (or other temperature in- 
dicators involving optical data) and near-IR properties, a 
regular behaviour is observed along the average sequence. 
This would not have been the case if sorting had been 
based on near-IR indices. For C-rich stars, both a tem- 
perature sequence and a sequence based on the C/O ratio 
are presented. For stars with thick dust envelopes (OH/IR 
stars and their carbon rich equivalent), the use of reddened 
versions of the above data is suggested. 

In studies of stellar populations, the average spectra 
will be used in connection with stellar evolution tracks 
or isochrones. The latter are expected to provide the dis- 
tribution of stars in the HR diagram, but also the on- 
set of LPV pulsation (i.e. when to use the average LPV 
spectra instead of static giant spectra), the transition be- 
tween O-rich and C-rich atmospheres, and the transition 
from an optically visible object to a dust-enshrouded far- 
IR source. As both the effective temperatures provided 
by stellar evolution tracks for the thermally pulsing AGB 
and the temperatures estimated for the average spectra 
are (independently) model-dependent and uncertain, we 
suggest caution when using the values indicated in this 
paper. The effects of reasonable changes in the tempera- 
ture scales have been illustrated by Langon et al. ( 1999| ). 
We recommend an a posteriori calibration of the relative 
temperature scales of the spectra and the evolutionary 
tracks, based on observed properties of intermediate age 
star clusters of known ages and metallicities. 

Using only one parameter, the effective temperature, 
to characterize LPV spectra is clearly a first order ap- 
proximation (even if this particular choice of a parameter 
is arguably the most appropriate one). Pulsation period 
and amplitude as well as the phase in the pulsation cy- 
cle are expected to produce systematic effects, although 
neither the FW2000 sample nor the currently available 
stellar models are sufficient to pin those down. The pe- 
riod and amplitude distributions in the LW2000 data used 
here are biased towards high values, when compared to the 
distribution of LPVs (including Miras and semi-regular 
variables) found in star counts in the Milky Way or the 
Magellanic Clouds. As only TP- AGB stars with relatively 



high main sequence masses will ever contribute much to 
the integrated light of stellar populations, the problem is 
not as severe there. In the practice of population synthesis 
calculations, the errors resulting from this period and am- 
plitude distribution cannot be disentangled from errors in 
the assumed onset instants of semi-regular and Mira-type 
pulsation. 
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